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TECHNICAL MEMORANDUM 1252 


. SOME STUDIES ON THE FLOW OF A GAS IN THE REGION OF TRANSITION 
THROUGH THE VELOCITY OF SOUND* 

By I. A. Kiebel 

The two-dimensional motion of an incompressible fluid about 
a closed contour with a definite velocity in magnitude and 
direction at infinity is considered. If, without changing the 
direction of the velocity at infinity, the magnitude is increased, 
the configuration of the streamlines remains unchanged and only 
the numbering of the stream function changes . There exists only 
one family of ourves that can serve as streamlines in the incom- 
pressible flow about a given contour (at a given angle of attack) j 
for example, the contour of an airplane wing. The case is quite 
different with a compressible fluid. For the components v x and 
Vy of the velocity along the x- and y-axes for incompressible 

flow 


T x _ a® a* 
v. ~ it * 3F 

T v dir 

"3E 

(voo is the magnitude of the velocity at infinity, qp the velocity 
potential, and i* the stream function) so that for i* there is 
obtained a simple Laplace equation. For a compressible gas moving 
adiabatically, there must be obtained 


^x = dcp = PO 

Vos p 5y 


*" Nekotorye Raboty o Teoheniakh Gaza v Oblasti Perekhoda 
Cherez Skorost Zvuka." Izvestia Akademii Nauk SSSR, No. 3, 1947, 
pp. 253-259. 
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i = dcp . p 0 34/ 


where 


P 


v Is "the ma g ni tude of the velocity, a* the critical, velocity, and 
k the ratio of specif io heats • The equation obtained for the stre am 
function vill now oontain as a parameter (^/a*) 2 , so that to one 

profile there will correspond an infinity of streamlines represent- 
ing this profile for different values of the velocity at infinity. 

The change in shape of the streamlines is particularly s har p 
when in the flow plane (still having a subsonic velocity at infinity) 
supersonic zones arise* A new specific difficulty now arises; namely, 
that the flow of a compressible fluid has two singularities as 
compared with the motion of an incompressible fluid. In the first 
place in the compressible fluid in contrast with the incompressible 
fluid. Infinitely large ve locities are impossible (the maximum 

possible velooity Is a* /\j( k+l)/(k-l)); in the second place In the 
supersonic flow of a gas, the stream tubes. In contrast to what 
holds for the incompressible fluid, expand with increasing velocity. 
This last circumstance leads to the fact that in the supersonic 
zone the streamlines will diverge and move away from the boundary, 
whereas In the subsonic zone the streamlines will converge on 
approaching the supersonic zone. It may be expected that for a 
given contour there will exist velocities at infinity for which 
both these laws cannot be satisfied. The finiteness of the velocity 
on the other hand leads to the fact that where the solution for 
the Incompressible flow gives Infinity velocities (for example, in 
the flow about an edge), the solution for the compressible flow 
either does not exist or the existing strealines do not form a 
sharp angle. Mathematically expressed, In the supersonic zone 
points and entire lines may appear on whioh the derivatives of the 
velocities become infinite . Such solutions of the equations of 
gas dy nami cs, although existing formally, have no physical meaning 
and oannot be realized. In these cases, the motion apparently so 
readjusts Itself that a line of strong discontinuity arises and 
the solution from the very beginning must be sought not in the form 
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of continuous Irrotational flow "but in the form of a motion in 
which there is a surface of discontinuity, the shape and position 
of which sire unknown initially and after passing through which 
the motion "becomes rotational . Tests show that such Jumps are, 
as a rule, actually found on the wing in the supersonic zone (the 
wing moves with subsonic velocity) . 

A complete theory of the transition through the velocity of 
sou n d has not, as yet, been developed. The problem of the transi- 
tion through the velooity of sound is one of the primary problems 
of aerodynamics. 

To N. E. Joukcwsky belongs the credit for the known modification 
of the method of Kirchhoff, where in the study of the motion of an 
incompressible fluid there are used in place of the coordinates x 
and y the magnitudes 0 (the angle of inclination of the velocity) 
and Z = In v # /v as the independent variables. 

Expressed in these variables, the equations for the stream 
function \l/ and the velocity potential cp (the motion is con- 
sidered irrotational) the following form is assumed: 


dcp cW 

33T “ - 3z 


dcp d\|f 

3z " & 

so that for the determination of \j/, the Laplace equation is 
obtained 

3z 2 

In the oase of a compressible fluid, a natural generalization 
of the variable Z is the function S of the velocity, employed 
by S. A. Christianovich in his fundamental work (reference 1) in 
the form 





1 - 


dv + constant 
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where 


M = 


v 

a 


a velocity of sound 


The equations of S . A. Christianovioh for the c ompre s s ible 
fluid may be written in the form 



( 1 ) 


where 0, as before, is the angle of inclination of the velocity. 


and 



For the supersonic region, it is convenient to introduce, as 
was done in another report of S. A. Christianovioh (reference 2), 
the functions X and a from the relations 

- 1 



1 


<7=1-0 from reference 2 . 
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and the equations for cp and \{/ then assume the form. 






aa ^ as 


The functions K, S, X > and c depend, in accordance with 
the Bernoulli law, only on the velocity, and 


aJk - t 


x S LzJk 

h 2 - tM 2 

h 2 - 1 7 


1 fl +, t / h - t\ h 

2 ln 1 - t \h + a) 


-2 = * + 1 
1 k - 1 


t =M - M 


Mfc? 2 


h 2 - 1 


c s h arc tan — - arc tan t 
h 


■ I^M 2 - 1 


where 
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The function A JiE through the intermediary of the parameter t 
depends on S; the function Vx~ through the intermediary of t depends 
on o. In figure 1, the values of S are in the direction of 
the positive horizontal axis and the values of Vk along the 
positive vertical axis; on the negative horizontal axis are the 
values of - a and along the negative vertical axis, -a/3c» As the 
velocities increase from 0 (M = 0) to a* (M = 1), S decreases 
from oo to 0. The value of Vk up to very small values of S 
will he very near unity and only at S = 0 does / \(k rapidly tend 
to zero. Further, when o increases from zero to its limiting 


value equal to 
to infinity (for 

inflection) . 


a » - — 
2k 


th e function will asymptotically tend 


A+ — 

k? k 


the curve has a point of 


When cp is eliminated in the subsonic region, the following 
equation is eliminated 


afy + afy + d in Afg ajjf 
Se 2 S 2 as 26 


and in the supersonic region 


cf£ + + d In A/X W . 0 

be 2 Z# 33 

On the basis of the behavior of the function a/k^ in equa- 
tions (l) VlC = constant may be set to very small values of S. 

The equation of Laplace is then obtained, which is the first approxi- 
mation of S. A. Christianovlch (references 1 and 3). In one of his 
later papers, Christianovich showed (reference 4) that for the super- 
sonic region, _t hat is, where equations (2) must be taken, starting 
from small a and further on, the function / \JX may be approximated 
in the form of pieces of ourves of the second order (first a para- 
bola and. then a hyperbola) . Equations (2) then become the equation 
of Darboux and may be simply integrated to the end. 

The analysis is, however, most difficult of all for the narrow 
region corresponding to the values of S and C near zero. This 
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region, is shown, to a magnified scale in. figure 2. This is the 
region of transition through the sound velocity, which was pre- 
viously discussed. It is first necessary to mention the investi- 
gation of F. I. Frankl (references 5, 6), who first for the region 
of transition separated out the principal term, in the equations 
for cp and and obtained for the principal term the equation 


of Tricomi. Frankl gives a profound analysis of the phenomenon 
of the transition through the velocity of sound and in particular 
investigated the condition at which in the Laval nozzle a trans- 
ition without the occurrence of a surface of discontinuity is 
possible. Recently his coworker, S. V. Falkovioh (references 7 
and 8), on the basis of Frankl 's work notes the fact that the 
function '\Jk about S = 0 can be devloped into a series of the 


form 


>\[k 


about 

3 rr 

IS + 


where bj_ depends on 



Falkovich then 


assumes near the transition line the following equation for ^ : 


dfy 1_ 34 1 

de 2 + 3s 2 + 3S cS 


(3) 


Introducing In place of 8 and S plane bipolar coordinates 
a and p from the relations 


a 


In 


32 + (S + S 0 ) 2 

le 2 + (s - s 0 ) 2 


p = arc tan 


2S 0 e 


6 2 +S 2 - S 0 2 


where Sq is the value of S corresponding to the velocity at the 

the point where 8=0 (the velocity at infinity), Fal k ovioh was 
able to construct a number of important particular solutions of 


Frankl introduced as independent variables 8 and 

’‘-(I*)* =(l 5 ) 5 - 
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equation (3). Oving to the work of Falkovich, it is possible 
to construct for the motion near the velocity of sound simple solu- 
tions having singularities of the type that arise in an Incompres- 
sible fluid in the problem of the flow with circulation. 


The occurrence of a surface of discontinuity in the supersonic 
region greatly complicates the solution of the problem. After 
passing through the surface of discontinuity, the motion becomes 
rotational - the velocity potential does not apply. For the 
stream function, the following complicated equation is obtained: 



where 


(4) 


1 d In* NO 
k - 1 d\|r 


where a funotion only of \J/ (constant in the irrotational 

case), is given by the value $ = — P ^ and the form of this 

funotion depends on the shape of the surface of discontinuity; 
P 0 /p and M have the previous value s 




3 . 


3 

As in the irrotational case, a* remains constant. 
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In studying the rotational motions, It would apparently 
he more convenient to proceed differently; name ly , to take as 
independent variables one of the coordinates, for example, x 
ana the str eam function \}/ . By making use of these independent 


variables and introducing as the required function 


pdx (the 


integration is taken for fixed \(/) , a number of particular exact 
solutions of the rotational problem (reference 9) are obtained. 

In particular, a solution is constructed with a transition through 
the velocity of so und - At the AN SSSR Institute of Mecha ni cs, 
investigations are being continued on the transition zone in both 
the rotational and irrotational oases . 


Translated by S. Beiss, 
National Advisory Committee 
for Aeronautics. 
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